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The scalar perturbations in inflationary models, based on a two-component diagonal non-linear sigma model, 
are considered. For inhomogeneities generated at an inflationary stage, the law of motion of the comoving 
curvature 1Z is obtained (without using the slow roll approximations). A formal expressions for the power 
spectrum and its spectral index are obtained, which are valid and after the slow-roll stage. As an example, 
an inflationary model with a massive scalar field and quadratic curvature corrections is studied numerically. 



^ ! 1 Introduction. 



A standard hypothesis of modern cosmology is the exis- 
tence of an accelerated expansion stage in the history of 
the very early Universe, when the second time derivative 
of the scale factor is positive. This expansion was termed 
inflationary. Presence of a long enough inflationary stage 
can resolve such well-known difficulties of the Big Bang 
theory, as the homogeneity, horizon and flatness problems 
P. Besides, the inflationary theory gives a root source of 
primordial inhomogeneities, giving rise to the large-scale 
structure of the Universe. In the simple models of chaotic 
inflation the existence only one scalar field termed inflaton 
is supposed. In a wide range of initial values of this field, 
the slow roll conditions is realized, and the scale factor 
grows quasi exponential. The primordial density inhomo- 
geneities arise from quantum vacuum fluctuations of the 
inflaton field. These quantum fluctuations, after horizon 
crossing, can be considered as classical stochastic pertur- 
bations PJ. At large scales, the conservation law of the co- 
moving curvature perturbation 1Z allows to connect inho- 
mogeneities generated on inflationary stage with primor- 
dial inhomogeneities in a radiation-dominated Universe. 
The slow-roll models predict initial inhomogeneities with 
an approximately scale-invariant (Harrison - Zel'dovich) 
spectrum in good agreement with the astrophysical data 
Pj. However, more realistic multi- field inflationary mod- 
els also call the attention. In such models, even for long- 
wavelength inhomogeneities, the conservation law of the 
quantity 1Z may not to hold (for example, strong ampli- 
fication of long-wavelength inhomogeneities is possible at 
the preheating stage). Many authors considered the evolu- 
tion TZ for a number of multi-field models both analytically 
and numerically. For double-field inflationary models, the 
equation of motion of the comoving curvature perturba- 
tion was written out in the pioneering papers 0,0. 

In this work, we consider inflationary models based on a 
non-linear sigma model with self-coupling potential. The 



general two-component diagonal sigma-model action has 
the form: 

s = |d 4 x/M{-^ + ^ 2 (^,xV; M ^ 



+ 2^te.x)X;MX iM -E%,x)j, (!) 

where the chiral metric components T 2 (tp, P 2 (<p, x) are 
some functions of chiral fields (sigma model components) 
if and x- The non-linear sigma models usually arise after a 
conformal transformation in generalized Einstein theories 
|B] . In Refs. , [Hj , inflationary models with an action of 
such view have received the name of chiral models. They 
are often call simply multi- field models |5JE], however, 
we do not use this term to stress the non-triviality of the 
scalar fields space. 

In two-field inflationary models, to study adiabatic and 
entropy perturbations, Gordon at al. jl lj constructed the 
quantities 5a and 6s . These variables are very convenient 
for describing the time evolution of the comoving curva- 
ture TZ, and for solving the perturbed field and Einstein 
equations in slow roll. The perturbation equations in the 
new variables are also suitable for numerical simulations 
I12| . This approach was generalized to the case of 
a non-linear sigma model with chiral metric components 
h n = 1, h u = 0, h 
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= P 26((pi) 



[T3| , and to the case of the 
general non-linear sigma model under slow roll conditions 
|10|. The goal of the presented work is to study the evolu- 
tion of comoving curvature perturbations and to obtain a 
power spectrum of initial inhomogeneities in inflationary 
models with the action (JIJ. 

2 Background. 

We consider the spatially flat Friedmann - Robertson - 
Walker (FRW) universe, described by the action (JIJ with 
the line clement 
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ds 2 = alt 2 - a 2 (t)S lj dx i dx j . 



(2) 
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The equations of motion for the homogeneous back- 
ground chiral fields <p, \ and scale factor a(t) are 
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where H = - is the Hubble rate. 

a __ __ 

Following the papers we define the new "adia- 

batic" held a such that 



a = Tip cos 9 + Px sin 9 , 



(7) 



where 



: —COS 9, 



Px 



^jT 2 p 2 + P 2 X 2 "~ ' y / T 2 L P 2 + P 2 x 2 
Using definitions 

U a = ^cos9 + ^sin9, 
[/ s = ^cos9-^sin9, 



= =sin9. (8) 



(9) 



from the chiral field equations (JjjJ and © we can write 
very useful relations: 



a + 3H& + U a = 0, 



UP T 

9 = —y H —in. 

& T P 



(10) 
(11) 



3 Perturbations. 



We investigate scalar linear perturbations about the spa- 
tially flat FRW background. The corresponding perturbed 
line element can in general be written as |141 115| 

ds 2 = (1 + 2(j>)dt 2 - 2aB. l dx i dt 

-a 2 [(1 - 2V>)4; + 2E tij ] dx l dx j . (12) 

The non-linear sigma model components can also be 
decomposed into a homogeneous background fields <p(t), 
x{t) and small inhomogeneities Sip, Sx- 



ip(x,t) = ip(t) +8<p(x,t), 



(13) 



Gauge transformations t = t + £°(t, x), x = x + £ ,l (i,x), 
where £°(t, x), £(t, x) are the arbitrary scalar functions, al- 
low one to simplify the metric tensor by a suitable choice 



of constant-time hypersurfaces and spatial coordinates in- 
side them. Here we operate without specified any partic- 
ular gauge conditions. 

Using perturbations of the energy-momentum tensor 
ST{f and those of the Einstein tensor 5G% , one can ob- 
tain the perturbed Einstein equations for inhomogeneities 
with the comoving wave number k: 



3#(#0 k + 4) + % 



V^k + H(a 2 E k - aP k ) 

K 2 

= — TrSpk, 



fa + 3PVk + #0k + (2H + 3P 2 )0 k 

K 2 



(14) 
(15) 

(16) 



0k- ^k+ P(aP k - a 2 Pk) + (oB k - a 2 P k )' = 0, (17) 
where the quantities Sp^, 5g k are 

S Pk = 1 a (p5(f> k + {U v + PP, v x 2 + TT^Sipk 
+P 2 xSXk + (U, x + PP, x X 2 +TT x ip 2 )S Xk 'CT 2 ^ (18) 

Sq k = -T 2 p6^ - P 2 xSxm- (19) 

In studying the evolution of perturbations, gauge- 
invariant quantities are very convenient. The following 
gauge-invariant quantities are often used: 



$ = 0+ (aB-a 2 E) , 



* = ip - H (aB - a 2 E 



(20) 



They were constructed by Bardeen ^H] and are widely 
used in the influential report |14j . The comoving curvature 
perturbation IZ^ can be expressed using the quantities $k 
and 4"k as follows |17| : 



ff (ff$ k + # k ) 
H 



(21) 



Eq. Q17[l can be rewritten in the simple form 

$ k = * k . (22) 

Taking a time derivative of Eq. 121(1 . using the back- 
ground field equations and Eqs. (fT^|) - JT7J|, we find the 
equation of motion of the comoving curvature perturba- 
tion 



K k = 2Hipx 



. .T 2 pU x -P 2 X U v fSp k <5 Xk 



{T 2 ip 2 + P 2 x 2 ) 2 \ ip X 

Hk 2 
+-■ — 
Ha 



2 * k . (23) 
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In the description of the evolution of the inhomo- 
geneities, we shall also use the new variables Ss and Sa 



Ss = P cos QS\ — T sin Q5<p, 
Sa = P sin Q5x + T cos QSip. 



(24) 



The quantities Sa and Ss were termed in Ref. "adi- 
abatic" and "entropy" perturbations. As indicated, by the 
name, perturbations with Ss = are only adiabatic. By 
definition, it is clear that Sa is associated with inhomo- 
geneities of the field a while Ss corresponds to inhomo- 
geneities of the entropy field s defined by 



ds = P cos Qd\ — T sin Od(p. 



(25) 



In some two-field inflationary models |18|. there is a very 
simple relationship between the generalized entropy per- 
turbations Ss after inflation and the specific entropy at the 
start of radiation dominated stage [111 IT§| . 

Using variables Sa and Ss, Eq. Il2.'i[l can be reduce to 



H h 2 TT 
H a a 1 



(26) 



which formally coincides with the corresponding equations 
of Refs. • For purely adiabatic perturbations 

(Sip/ip = §x/X)i we have Ss = 0, and at large scales 
[k/a <C H), a conservation law is valid for comoving cur- 
vature perturbation TZ^. 

The perturbed chiral field equations also can be pre- 
sented in the terms of the new variables. Directly per- 
turbing the chiral field equations [201 for sigma model QJ, 
we find: 



^ k +(3tf+2^W k +2fe- xW 
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After some calculations, the following equations for the 
inhomogeneities Sa^ and <5s k can be derived: 



I k 2 U ■ UP 
Sa^ + 3HSa^+ \— + -LQ + U aa -^—. ^sinO 
\ cl a a P z 



^—cosQ\Sa k + 2[—Ss k ) -2- — -Ss v 



a 



a a 

.2 i 



fc Pw 

= -2P CT k + & ( k + 3^ k + k 2 E k - " 



a 



Ssk + 3HSsk 



P T 
U 



PT 



(P lV> cos 6 + T x sin 6) 



P.L 



T, 



p2 



-r— | C0S6 
P A 

sine 



T 
T 



r 



sine 



P 



a 

PT 



T T 
— sin 6 # cos e 



5> 
P 

p 



• cos e 



.17. 



where 
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(29) 



Here we use the gauge-invariant quantity e m 
3HSqk the comoving density perturbation: 



(30) 
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(32) 
Sp-k - 



,(a 3 TT x ip) (T 2 )^ x , (P 2 ), yx ,o 



a 3 T 2 

fc 2 t/, 
^ + ^T 2 



2T 2 



2T 2 



X \5xk 



^<Vk + ^Xk + 2^f< 



rp2 



<p k + 3^ k 



fc 2 Pk , 



( 2 (^gxX) (^x X ,, 2 (T 2 ), xx 



a 3 P 2 



2P 2 



a 3p2 



2P 2 



2P 2 
2P 2 



<^ 2 ^Xi* 



(27) 



<5(j k + 2—^-Ssk — (T0k — -r<5crk 
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Since from the Eqs. (|14|) and l|15|) we have 



k 2 K 

or 2 



(33) 



(34) 



at large scales (fc /a <C H) the right hand side of Eq. (|3*U)l 
becomes negligible. As follows from Eq. H3U|) . the large- 
scale entropy perturbations are decoupled from adiabatic 
and metric perturbations, which generalizes the conclusion 
of Ref. on inflationary models with the action QJ . 

Eq. Q29J) can be rewritten as an equation for Sasaki- 
Mukhanov gauge-invariant variable 



Sa + —ip. 



(35) 
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Using this quantity, one can obtain: 

Q k + 3#Q k + (^ + %e-%4sine 



U t9 T n 2 fa 3 a 2 

— — — K cos B - 

a T a 3 \ H 



Eqs. (HU and ® in specific case T 2 = 1, P 2 = P 2 {tp) 
coincide with the corresponding equations of Ref. |13| . 
and in the slow roll conditions reduce to equations of Ref. 

D33- 

4 Application to slow roll infla- 
tionary models. 

Let us consider an application of the above expressions to 
slow roll inflationary models. Quantum vacuum fluctua- 
tions of chiral fields become classical inhomogeneities after 
they leave the horizon. The Fourier component of chiral 
fields perturbations at horizon crossing are given by |S] 



<fy>k(*k) = 



1 ff(t k ) 
T(i k ) V^k 3 

1 H(t k ) 
P(tk) 



e<p(k), 
e x (k). 



(37) 
(38) 



Here ik is the instant of horizon crossing (a(tk)-ff(ik) = k), 
and e^, e x are real Gaussian stochastic variables with the 
following properties: 



Mk)e x (k')) = <W5 (3) (k-k'), 
Mk)) = (e x (k)> = 0. 



(39) 
(40) 



From l|24|l and (|37|l one can derive the following expres- 
sions for dak and Ss^: 

<M*k) = ^=e CT (k), Ss k (t k ) = ^L(k), (41) 



/2fc 3 



'2k 3 



where e CT , e s are independent Gaussian random quantities 
with zero average and unit dispersion. 

Being restricted to the most important non-decreasing 
modes of long-wavelength inhomogeneities, as well as in 
case of single-field inflationary models, it is possible to 
write the inequality: 



*k « #*k- 



(42) 



Another standard assumption is that, for non- 
decreasing modes in the slow-roll regime \5ip/5ip\ <C 3H, 
\S'x/S'x\ < 3iJ and therefore Q3|, EH 



\Sa\ < 3H\6a\, 



\6s\ < 3H\5s\ 



(43) 



This allows us to disregard Sa and da in the perturbed 
field equations J2HJl, Thus, Eq. (J20J may be written 

as 

S's k = A(t)Ss k , (44) 

where 

A ® = ~^(^ (P V cos 9 + T x sin 6) 
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Taking into account the inequality l|42|) . equations (|15|l 
and l|21|) can be simplified: 



2i? 



(T(5cr k , 



H 2 

n k = — p -^> k . 

H 



(46) 
(47) 



Integrating Eq. I|26(l . one can obtain the following ex- 
pression for the long-wavelength comoving curvature per- 
turbation IZk : 



ftk(t) = ^k(ik) 



g(*k) 
V2J 3 



J(t,i k )e s (k), 



where 
J(i, tk) 



2H£ 1 _Ss_ dt _ 
o 2 Ss(t k ) Jt 



— r^—aUtjdt. 
a z 



(48) 



(49) 



Here the function atk(t) is a solution of Eq. (|44() . nor- 
malized by the condition ak(£k) = 1- Using first of re- 
lations l|41l) . background equation Q and also equations 
l|46|l . (|47|l one can therefore write: 



#(*k) /^(*k 



2fc 3 \<K*k) 



e CT (k)- J(i,i k )e s (k) 



(50) 



The power spectrum |15) can now be written in the fol- 
lowing form: 



H 2 (t k ) fH 2 (t k 



4tt 2 



a 2 (t k ) 



(51) 



This formula is valid also after the slow-roll stage if, 
for finding the function a k (i) a long- wavelength limit of 
Eq. ll.'iUII is used instead of Eq. I|44|l . If, from certain time, 
the entropy perturbation Ss is so small, that the comoving 
curvature perturbation 1Z remains actually constant, the 
expression obtained above describe the adiabatic power 
spectrum. 

A commonly used characteristic of a power spectrum 
is its spectral index n s ^BJ. Using l|51() and background 
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Figure 2: Numerical simulations of the spectral index af- 
Figure 1: The time dependence of V(t) /V(t k ) with ter inflation. Here it is assumed that the perturbations 



M/m = 5, <p(t k ) = 14/k, X (tk) = 4//s. 



equations one can conclude, that the spectral index for 
inflationary models under consideration is given by: 



n l _ d\nV n {t) 
d\nk 

J2 , 2g(t k )q(t k ) 

K + <T 3 (« k ) 



1 d\nV n {t) =2 _H(t k ) 



H 



dt k 

' 2U s {t k ) 
v * 2 (*k) ' 



H 2 (t k ) 



A(t k ) 
H(t k ) 



H 2 (t k ) 



J 2 



(52) 



The spectral index is time-dependent due the time- 
dependence of the quantity J(t, t k ). 

Eqs. jnU), J2E1), (EDI and J3fJ) are the main results of this 
work. These equations are very well suitable for numer- 
ical calculations of perturbation evolution in inflationary 
models based on non-linear sigma models. 

4.1 Numerical example. 

As an example, we shall describe an inflationary model 
with a single massive scalar field and additional higher 
order curvature terms due to curvature squared. The cor- 
responding action looks like: 



S 



d 4 x 




2 2 

m if 



(53) 



Making the conformal transformation g^ — fl g%j with 



n 2 



gjp-, one gets the action in the Einstein frame: 
R 1 .„ e-^ KX 



S = I d xy/\g\<— + -x- tl x 



3M 2 



T K V 



\e- 2 \fl«*m 2 v 2 \, (54) 



leave the horizon 60 e-folds before the end of inflation; the 
inflation ends when H/3H 2 = 1. 



where the new field \ = In (f - R 2 /3M 2 

introduced. Note that when \ then U(ip,x) ~ 

M 2 x 2 /2+m 2 ip 2 /2. Some inflationary models more general 
than were considered earlier in the preheating context 
in (211122) (These papers took into account not only higher 
curvature terms, but also a possible non-minimal coupling 
& 2 R). 

We are interesting here in the case when the chiral fields 
if and x are comparable in magnitude and are both in a 
slow-roll regime at the beginning of the inflationary stage. 
In the analysis of the background and inhomogeneities 
evolution, we shall be restricted to the parameter area 
M > m. Since the field x 1S heavier, it breaks up faster 
than field ip, and soon after the inflation influence of the 
field x on the scale factor dynamic and evolution of inho- 
mogeneities is vanishingly small. Besides, after inflation, 
f2 2 = 1, therefore actions in the original frame l|53|l and 
in the Einstein frame l|54(l almost do not differ. Accord- 
ingly, the inhomogeneities calculated in both frames after 
inflation asymptotically coincide. It is clear that the inho- 
mogeneities in this model after the inflationary stage are 
to a great extent adiabatic. 

The background and perturbed equations for the sigma 
model (I54|l can be solved by numerical simulations. For 
finding the long-wavelength comoving curvature perturba- 
tion 1Z generated on inflationary stage, it is very conve- 
nient to use Eqs. H26J) and 1(30(1 . Eq. 1(30(1 is much easier 
to solve than the initial set of perturbed Einstein and field 
equations. It can be made by using standard packages of 
numerical computations. The results of the calculations 
are presented in Figs. 1 and 2. These pictures show the 
time dependence V(t)/P(t k ) for a specific choice of pa- 
rameters and the dependence of the spectral index on the 
model parameters. 
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